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C/^ ' Abstract 

We study the dynamics of a 'composite' or 'dressed' quark in strongly- 
coupled large- A'^c = 4 super- Yang- Mills, making use of the AdS/CFT cor- 
respondence. We show that the standard string dynamics nicely captures the 
^ • physics of the quark and its surrounding non-Abelian field configuration, mak- 

1^ , ing it possible to derive a relativistic equation of motion that incorporates the 

effects of radiation damping. From this equation one can deduce a non-standard 
^^^^ dispersion relation for the composite quark, as well as a Lorentz covariant for- 

^ . mula for its rate of radiation. We explore the consequences of the equation in 

' a few simple examples. 

a^ ■ 
in 

^ , 1 Introduction and Summary 

o . 

^ . When a charge radiates, energy conservation dictates that it must be subjected to a 

reactive force originating from its self-field, that tends to damp its motion. In the 
context of classical electrodynamics, the study of this damping or radiation reaction 
^ ■ force began over a century ago [H El [3l H] , and continues to this day [5], [6]. Reviews 

■ and additional references on the subject may be found in [3, El El fTO] . 

In a non-relativistic approximation, the dynamics of an electron that is modeled 
as a vanishingly small spherically symmetric charge distribution is controlled by the 
Abraham-Lorentz equation [H [2] 

m (^x — teX 

where ' = d/dt and tg = 2e^/3mc^ is a timescale set by the classical electron radius. 
In this equation, the damping force (the second term in the left-hand side) is seen to 
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be proportional to the jerk J = a = v. The search for a Lorentz-covariant version of 
([T]) led to the (Abraham-) Lorentz-Dirac equation [1], 



m ( F — t 



, (2) 



with ° = d/dr, r the proper time (defined such that F = — c^) and J-"^ = 
^{F -v/c, F) the four- force. The second term within the square brackets (proportional 
to the square of the proper acceleration) is the negative of the rate at which energy and 
momentum is carried away from the charge by radiation, according to the covariant 
Lienard(-Larmor-Heaviside-Abraham) formula. So, strictly speaking, it is only this 
term that can properly be called radiation reaction. The first term within the square 
brackets, usually called the Schott term, and whose spatial part yields the damping 
force of ([1]) in the non-relativistic limit, is known to arise from the effect of the charge's 
'near' or 'bound' (as opposed to radiation) field pTl 17]. 

The appearance of third-order terms in ([T]) and ([2]) leads to unphysical behavior, 
including pre-accelerating and self-accelerating (or 'runaway') solutions. These defi- 
ciencies are known to originate from the assumption that the charge is pointlike0 For 
a charge distribution of small but finite size /, the above equations can be shown to be 
truncations of expressions that involve an infinite number of derivatives {Id/dt)"- (and 
generally include terms that are non- linear in these derivatives), but are physically 
sound as long as / > ctg and [3, [3 [9] . 

Of course, one should keep in mind that the unphysical behavior implied by ([1]) 
and ([2]) would be visible only for time and distance scales smaller than the Compton 
wavelength Ac = h/m of the charge, and thus lies outside of the actual range of 
validity of classical electrodynamics. How the preceding story generalizes to the case 
of fully quantum electrodynamics (QED) has been studied from different angles in 
[T2| [T3l [HI [15] and references therein. In particular, in [12] it was shown that, for 
a pointlike non-relativistic electron, QED leads to an equation of motion with an 
infinite number of higher derivatives, implying that the electron acquires an effective 
size / = \c due to its surrounding cloud of virtual particles. 

Going further to non-Abelian gauge theories is a serious challenge!^ paper to show 
that the AdS/CFT correspondence [HI [19], [20] allows us to examine this question 
rather easily in quantum strongly-coupled non-Abelian gauge theories. The essence 
of the matter is that in the context of this duality the quark corresponds to the tip 
of a string, whose body codifies the profile of the non-Abelian (near and radiation) 
fields sourced by the quark. In other words, the quark has a tail, and it is this tail 
that is responsible for the damping force. Indeed, this mechanism has already been 
seen at work in the recent computations of the drag force exerted on the quark by a 



^This assumption leads to the further comphcation of an infinite electromagnetic self-energy, 
which has already been been absorbed within the renormalized mass m shown in (jT]) and An 
alternative approach which circumvents this divergence has been proposed recently in [S] . 

^Scc respectively }16| and [17j for work on radiation within classical and (weakly-coupled) quan- 
tum Yang-Mills theory. 



2 



thermal plasma, which is described in dual language in terms of a string living on a 
black hole geometry [2T],[22]. Our analysis makes it clear that, irrespective of whether 
a spacetime black hole is present or not, the body of the string plays the role of an 
energy sink, as befits its identification as the embodiment of the gluonic degrees of 
freedomjfl 

We expect this basic story to apply generally to all examples of the gauge/string 
duality, including cases with finite temperature or chemical potentials, but for sim- 
plicity we will concentrate on the case of quark motion in the vacuum of A/" = 4 
super- Yang-Mills (SYM), where, building upon previous work [301 [25], we can achieve 
full analytic control. As has been remarked several times in the past, it is interesting 
that even in this non-confining theory the gluonic field configuration can be encoded 
in a 'QCD' string, albeit one that lives in a curved higher dimensional spacetime. 
Our results are a direct consequence of this amazing fact. 

The paper is organized as follows. In Section [2], we explain how the standard 
string dynamics is mapped by the AdS/CFT correspondence onto the dynamics of 
an extended radiating particle. We begin by setting up our problem in Section 12. ![ 
reviewing the relevant context and emphasizing some key features — most notably, the 
fact that the quark with finite mass that the correspondence puts at our disposal is 
automatically 'dressed' or 'composite', as discussed around Eq. Qj. We then attack 
the problem in Section 12. 2[ where we derive an equation of motion for this quark, 
Eq. (128|) . which constitutes our main result]^ As explained at length in the paragraphs 
that follow it, this equation is a nonlinear generalization of ([2D, which incorporates 
the effects of radiation damping, but has no pre-accelerating or self-ac celerating 
solutions. From this equation one can read off a non-standard dispersion relation 
for the quark, Eq. flMj) . as well as a Lorentz covariant formula for its radiation rate, 
Eq. (l35l) . An interesting novel feature in these expressions is their dependence on the 
external force exerted on the quark, which is a reflection of its extended, and hence 
deformable, nature. We close Section [2] by commenting on our failure to rewrite (1281) 
in terms of an action principle. 

In Section [3l we explore some of the physics implied by (!28l) , specializing to a 
few simple examples. For the case of one- dimensional motion, we show in Section 13.11 
that even though, as expected on physical grounds, zero external force implies zero 
acceleration, the converse is not true: the quark will not accelerate when subjected 
to an external force that takes the specific form (1411) (which is identically zero only 
when the parameter to ±oo). More generally, for each given quark trajectory there 
is a one-parameter family of possible external forces. This again is a manifestation 
of the fact that, because of the extended character of the quark, the energy supplied 
to it can not only increase its velocity, but also modify its associated gluonic field 
profile. We end the paper by studying the nonrelativistic limit of fl2S]) in Section 13. 2^ 

■^On the other hand, energy loss via the strmg does turn out to be closely associated with the ap- 
pearance of a worldsheet horizon, as noticed initially in [23(l24j at finite temperature and emphasized 
in [25] for the zero temperature case. This association has been further studied in [26 l [27 | [28 l [29] . 
brief report of this derivation was given in the recent letter |31j . 



3 



where the hnearized form of the expressions allows us to make direct contact with 
the energy analysis of ([501 ESj) and to easily write down an action principle. 

All in all, then, we have in (1281) a physically sensible and interesting description 
of the dynamics of a composite quark in A/" = 4 SYM. This result serves, on the one 
hand, to illustrate the power of the AdS / CFT correspondence, and on the other, to 
shed some light on the largely uncharted terrain of radiation in strongly-coupled non- 
Abelian gauge theories. It would be interesting to extend this analysis in a number 
of directions. In particular, it seems worthwhile to explore the manner in which the 
split between intrinsic and radiated energy (and momentum) of the quark achieved 
in [30l [25] and the present paper, via examination of the string worldsheet, manifests 
itself in the gluonic field profile, by directly computing the expectation value of the 
energy-momentum tensor or similar local operators [32]. It is also natural to try to 
carry over some of the present methods to the finite temperature context [53] , where 
one sh ould be able to make contact with previous AdS/CFT analyses of energy 
loss in a thermal plasma (a rather large body of work detonated by the seminal 
works [211 [22l [Ml [35]), including the interesting recent studies of Brownian motion 
[5Sl[57l[5B]. 

2 From Strings to Quarks 
2.1 Basic setup 

It is by now well-known that strongly-coupled TV = 4 SU{Nc) SYM with coupling qym 
is dual to Type IIB string theory on a background that asymptotically approaches 
the AdSs x geometrjQ 

ds^ = GuNdx^Ux^ = — (-df + dx^ + dz^) + R^dn^ , (3) 

2 

(with a constant dilaton and N^. units of Ramond-Ramond five-form flux through the 
five-sphere), where Is denotes the string length [18]. The radial direction z is mapped 
holographically into a variable length scale in the gauge theory [39]. The directions 
x^ = (t, x) are parallel to the AdS boundary 2; = and are directly identified with the 
gauge theory directions. The state of IIB string theory described by the unperturbed 
metric ([3]) corresponds to the vacuum of the A/" = 4 SYM theory, and the closed string 
sector describing (small or large) fluctuations on top of it fully captures the gluonic 
(-|- adjoint scalar and fermionic) physics. 

From the gauge theory perspective, the introduction of an open string sector 
associated with a stack of Nf D7-branes in the geometry ([3]) is equivalent to the 
addition of Nf hypermultiplets in the fundamental representation of the SU{Nc) 

^From this point on we work in natural units c = 1 = 
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gauge group, breaking the supersymmetry down to A/" = 2 ^0]. These are the degrees 
of freedom that we refer to as 'quarks,' even though they include both spin 1/2 and 
spin fields. For Nf -C N^, the backreaction of the D7-branes on the geometry can 
be sensibly neglected; in the field theory this corresponds to working in a 'quenched' 
approximation which disregards quark loops (as well as the positive beta function 
they would generate). The D7-branes cover the four gauge theory directions, and 
extend along the radial AdS direction up from the boundary at z = to a position 
z = z„i where they 'end' (meaning that the C that they are wrapped on shrinks 
down to zero size), which is inversely proportional to the quark mass, 

^ (4) 



27rm 



An isolated quark is dual to an open string that extends radially from the D7- 
branes to the AdS horizon at 2; ^ oo. The string dynamics follows as usual from the 
Nambu-Goto action 

1 



-S'ng = / rfV V^det^ = y d^a C^q , (5) 

where Qah = daX^dbX'^GMNi^) {o,,b = 0,1) denotes the induced metric on the 
worldsheet. In our work the entire string will be taken (consistently with the corre- 
sponding equations of motion) to lie at the 'North Pole' of the (the point where 
the S'^ C that the D7-branes are wrapped on collapses to zero size), so the angular 
components of the metric (which are associated with the orientation of the gauge 
theory fields in the internal SU (4) symmetry group) will not play any role, and the 
lower endpoint of the string will necessarily lie at 2 = 2;^- 

We can exert an external force F on the string endpoint by turning on an electric 
field Foi = Fi on the D7-branes. This amounts to adding to the Nambu-Goto action 
the usual minimal coupling 

Sf = J rfrA^(X(r,^„))9.X'^(r,2™) , 

or, in terms of the quark worldline, 

5f= / drA^{x{T))x''{T) . (6) 



Notice that the string is being described (as is customary) in first-quantized lan- 
guage, and, as long as it is sufficiently heavy, we are allowed to treat it semiclassically. 
In gauge theory language, then, we are coupling a first-quantized quark to the gluonic 
(+ other SYM) field(s), and then carrying out the full path integral over the strongly- 
coupled field(s) (the result of which is codified by the AdS spacetime), but treating 
the path integral over the quark trajectory x^(r) in a saddle-point approximation. 
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Variation of the string action S'ng + 'S'f implies the standard Nambu-Goto equation 
of motion for all interior points of the string, plus the standard boundary condition 



where 



^ dCj<sG 




^{drX ■ d^xy - {drxni + {d^xy 



(7) 



is the worldsheet (Noether) current associated with spacetime momentum, and we 
haver recognized JF^ = —F^^drX" = (— 7-F ■ v,^F) as the Lorentz four-force. 

For the interpretation of our results it will be crucial to keep in mind that the 
quark described by this string is not 'bare' but 'composite' or 'dressed'. This can be 
seen most clearly by working out the expectation value of the gluonic field surrounding 
a static quark located at the origir|§ [13], 



(TrF2(x)) 



167r2|x|4 



T j_ 5 f 2nm\x\ 



5/2 



(9) 



For m — ^ cxD (zm — ^ 0), this is just the Coulombic field expected (by conformal 
invariance) for a pointlike charge. For finite m the profile is still Coulombic far away 
from the origin but in fact becomes non-singular at the location of the quark. 



4^ 



2 (TrF2(x)) 

YM 



1287r2 
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35 / 27rm|x| 



\x\ 



+ ... 



for Ixl < 



27rm 



As seen in these equations, the characteristic thickness of this non-Abelian charge 
distribution is precisely the length scale Zm defined in (jlj). This is then the size of the 
gluonic cloud that surrounds the quark, or in other words, the analog of the Compton 
wavelength for our non-Abelian source. 

It is interesting to note that the string can be alternatively viewed as a Born- 
Infeld string, i.e., a soliton of the gauge and scalar fields on the D7-brane [Hj. Since 
the small fluctuations of these fields (corresponding to microscopic open strings) are 
known to be dual to mesons, the composite quark itself can be thought of as a soliton 
constructed by aligning a large number of mesons [15]. The cloud surrounding our 
quark is then best thought of as 'mesonic' rather than 'gluonic'. Mesons are indeed 
known to be the lightest states in the spectrum of the strongly-coupled gauge theory, 
with masses of order mmcs = 1 /-^m = 2-7rm / <C m [l6] , and form factors with size 
set by Zm [IZ]. 



^More precisely, the operator in the left-hand side of ([9]) is the dual of the dilaton field, and 
includes not only the standard Yang-Mills term but also scalar and fermion contributions that can 
be found in [42] , and which we suppress for notational simplicity. 
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So, to summarize, Zm can properly be called the quark Compton wavelength in- 
sofar as it gives the size of the cloud of virtual particles surrounding the quark, but 
one should bear in mind that it is given not by 1/m but by l/mmes, and in this sense 
it could also be referred to as the meson Compton wavelength. 

2.2 Equation of motion for the quark 

The first analysis of an accelerating quark via the AdS/CFT correspondence was 
carried out in [H], which used tools developed in [19] to study the dilatonic waves 
given off by small fluctuations on a radial string in AdSs, and infer from them the 
profile of the gluonic field (TrF^(x)) in the presence of a quark undergoing small 
oscillations. The results of |1H] painted an interesting picture of the propagation of 
nonlinear waves in A/" = 4 SYM, but did not allow a definite identification of waves 
with the l/\x\ falloff associated with radiation. (More recently, this falloff has been 
successfully detected in the same setup through a calculation of the energy- momentum 
tensor {T^^) [50].) 

In [IH] it was noted that for Zm — > and in the linearized approximation, the string 
action ([5]) correctly implies the expected action for an ordinary non-relativistic particle 
of mass m — >■ oo. In the present section we will obtain the relativistic generalization 
of this result retaining the full non-linear structure of the Nambu-Goto string, and 
then further extend the analysis to the case of finite m. 

We will take as our starting point the results obtained in a remarkable paper 
by Mikhailov [3D], which we now briefly review (a more detailed explanation can be 
found in [23]). This author considered an infinitely massive quark, and was able to 
solve the equation of motion for the dual string on AdSs, for an arbitrary timelike 
trajectory of the string endpoint. In terms of the coordinates used in ([3]), his solution 
is 

X'^(r,^) = z^^^-^ + x'^(r) , (10) 
dr 

with a;'^(r) the worldline of the string endpoint at the AdS boundary — or, equiv- 
alently, the worldline of the dual, infinitely massive, quark — parametrized by its 
proper time r. 

Combining and flTU]) . the induced metric on the worldsheet is found to be 

9tt = ^(^^3-2 - 1), 9zz = 0, gzr = 7' (11) 

Z"^ z^ 

implying in particular that the constant-r lines are null, a fact that plays an important 
role in Mikhailov's construction. The solution (fTOl) is 'retarded', in the sense that the 
behavior at time t = X^{t,z) of the string segment located at radial position z 
is completely determined by the behavior of the string endpoint at an earlier time 
tj-etit, z) obtained by projecting back toward the boundary along the null line at fixed 
r. An analogous 'advanced' solution built upon the same endpoint /quark trajectory 
can be obtained by reversing the sign of the first term in the right-hand side of ffTOj) . 
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In gauge theory language, this choice of sign corresponds to the choice between a 
purely outgoing or purely ingoing boundary condition for the waves in the gluonic 
field at spatial infinity. Both on the string and the gauge theory sides, more general 
configurations should of course exist, but obtaining them explicitly is difficult due to 
the highly non-linear character of the system. Henceforth we will focus solely on the 
retarded solutions, which are the ones that capture the physics of present interest, 
with influences propagating outward from the quark to infinity. 

From the /j, = component of flTUl) . parametrizing the quark worldline by x^{t) 
instead of r, and using dr = — v'^dx^, where v = dx/dx^, the relation that defines 
the retarded time follows as 

^ + tret , (12) 



V 



2 



where the endpoint velocity v is meant to be evaluated at tret- In these same terms, 
the spatial components of ( ITOl) can be formulated as 

-> V 

X{t, Z) = Z + X{tret) = {t- tvet)v + f (tret) • (13) 

Using fll2p and (IT^ . Mikhailov was able to rewrite the total string energy in the form 

\/a /"* ,. d^ - \v X df 



2vr (1 - v^) 

where of course a = dv/dx^ . The first term codifies the accumulated energy lost by 
the quark over all times prior to t, and is surprisingly seen to have precisely the same 
form as the standard Lienard formula from classical electrodynamicsjll The second 
term in the above equation arises from a total derivative on the string worldsheet, 
and gives the expected Lorentz-covariant expression for the energy intrinsic to the 
quark [25] . 

/ 1 1^ "'"=° 



2tt \^i^rwz 

For the spatial momentum, [30| 125] similarly find 



7m . (15) 



a/a /"* ,. - f-u X al^ 



P{t) = IT / dirct—-^—-^v + pMt)) , (16) 



with 



'-ymv . (17) 



27r \^iZr^z, 

We see then that, in spite of the non-linear nature of the system, Mikhailov's pro- 
cedure leads to a clean separation between the tip and the tail of the string, i.e.. 



^Possible experimental implications of this result have been explored in |51j 
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between the quark (including its near field) and its gluonic radiation field. We will 
now exploit this separation to study in more detail the dynamics of the quark. 

Our initial observation is that, when we regard the Nambu-Goto action as a 
functional of the quark trajectory by plugging ffTTl) back into ([5]) + ([6]), we can 
explicitly carry out the integral over z to obtain 



which is evidently the standard action for a pointlike externally forced relativistic 
particle (with mass m — *■ oo). Notice that the associated equation of motion does 
not include a damping force, which is just as one would expect for an infinitely 
massive charge, because the coefficient te c(. 1/m oi the damping terms in ([T]) and 
approaches zero as m cxd. 

Let us now consider the more interesting case of a quark with finite mass, Zm > 0, 
where, as we emphasized in the previous subsection, our non-Abelian source is no 
longer pointlike but has size Zm- In this case the string endpoint is at 2; = z^,, and 
we must again require it to follow the given quark trajectory, x^{t). As before, this 
condition by itself does not pick out a unique string embedding. Just like we discussed 
for the infinitely massive case below (fTTl) . we additionally require the solution to be 
'retarded' or 'purely outgoing', in order to focus on the gluonic field causally set up by 
the quark. As in [25], we can inherit this structure by truncating a suitably selected 
retarded Mikhailov solution. The embeddings of interest to us can thus be regarded as 
the z> Zm portions of the solutions flTU]) . which are parametrized by data at the AdS 
boundary 2; = ol^ Henceforth we will use tildes to label these (now merely auxiliary) 
data, and distinguish them from the actual physical quantities (velocity, proper time, 
etc.) associated with the endpoint /quark aX z = Zm, which will be denoted without 
tildes. 

In this notation, (ITU]) reads 



Repeated differentiation of this equation with respect to f and evaluation at z = Zm 
(where we can read off the quark trajectory x^(f) = X^(f, Zm)) leads to the recursive 

*That this direct truncation indeed retains the retarded structure of the solutions is manifestly 
confirmed in our final rewriting (j29[) . where information is seen to propagate upward along the body 
of the string, i.e., from the UV to the IR of the gauge theory. 






(19) 
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relations 



IF 
~d^ 



z. 



d'^x^ dx^ 

1 — 

dr ' 

d'^x>' 



z. 



df^ 
d^x^^ 



dr 



•f2 



Adding these equations respectively multiplied by {—z^ 



(20) 



\n-l 



we can deduce that 



dx^^ dx^ 



dr dr 

and, upon further differentiation, 

d^x^" _ d^xf" 
df"^ df"^ 



d'^x^' 
'~d^ 



+ z' 



d^x^ 
'~d^ 



2 



(21) 



(22) 



This last expression already takes us halfway towards the equation we are after, but 
we still need to find a relation between df and the endpoint /quark proper time dr, 
and similarly rewrite d'^x^ /df"^ in terms of quantities at the actual string boundary 
z = Zm instead of the auxiliary data at 2 = 0. 
The first task is easy: from f|T9l) it follows that 



dX^ = dz 



dx^^ 
df 



dr 



d'^x'^ dx^ 
+ 



dr 



^2 



dr 



which evaluated at fixed z 



implies 



, „ ,^ , d^x^ dx^^ 
dx^ = dr I Zm, + 



dr 



^2 



dr 



and therefore 



dr = —dx^dx^ = df 



-2 



\-zL\ 



( d'^x 



\df^ 



(23) 



To arrive at this last equation, we have made use of the fact that r is by definition 
the proper time for the auxiliary worldline at z = 0, so (dx/df)"^ = — 1 and (dx/df) ■ 
(d^x/df^) = 0. 

What remains then is to express d'^x'^/df'^ as a function of quark data. For this 
we note first that, upon substituting the solution (fT9l) . the momentum current ([8]) 
(with appropriate tildes) evaluated at z = Zm simplifies to 



2n ~ 
W 



1 d^x., 



dr 



^2 



d x\ dx^ 



dr'^ J dr 



(24) 
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To avoid possible confusion, we should note that the tilde in the left-hand side does 
not indicate evaluation at z = (as all other tildes do), but the fact that this current 
is defined as charge (momentum) flow per unit f. The corresponding flow per unit 
r is clearly jusll^ 11^ = (9f/9r)n^, and it is this object which according to ([7]) must 
equal the external force JF^. Using this, fl23l) and the first equation of (j20|) in 
one can deduce (after some straightforward algebra) that 



^ _ 3 ^2 ^X11 



(25) 



where we have used the abbreviation = {2ti / \/X)T^. Since T'^dx^/dr = (which 
is merely the statement that no work is done on the quark in its instantaneous rest 



frame), this implies that [d'^x/dr'^Y 

df- 



j^^, which allows (123|) to be simplified into 



dr 



Using 
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and (1261) . we can finally rewrite (122]) in the form 



(26) 



dx^ 



d 



d 

dr 



dr 



dx^ 
dr 

dx^ 
dr 



(27) 



+ 



This equation of motion for the quark contains an infinite number of derivatives of 
x^, precisely as one would expect for an extended color charge distribution, based on 
the classical or quantum electrodynamic analogs [TJ El [12] mentioned in the Intro- 
duction. Notice that to arrive at this result we have made no assumption about the 
profile of the charge distribution. The dual string dynamics automatically incorpo- 
rates the physics of this profile, which is codified by the slope s = dzX{zm,t) [25]. It 
would be interesting to explore this connection in more detail through a calculation 
of {Tt F'^{x)) and similar observables for an accelerating quark in vacuum [IHl EDI [32] • 
A more manageable form of the equation of motion can be obtained by going back 
to the second equation in ([201) . 



d'^x" 



+ z„ 



d^x^ 



Through ([25]). ([26]) and 



df"^ dr^ dr^ 
this can be reexpressed as 




dr l-Kva 



' dr 



2nm 



A 



•772 
4-' 



f28) 



'The transformation rules for 11° under more general reparametrizations can be found in, e.g. 
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which is the equation we advertised in the Introduction. Notice that it involves only 
the four- velocity and four-acceleration of the quark, so, in going from (1271) to (125]) . we 
have traded an infinite number of higher derivatives for a somewhat more complicated 

dependence. This is to some extent analogous to the possibility of trading ([1]), ([2]) 
or its non-pointlike generalizations for an integro-differential (nonlocal in the force) 
equation with derivatives of only up to second order [71 [8] . What is different is 
that our end result, equation fl2Sl) . does not involve any nonlocality. (It is also highly 
nonlinear, because it includes effects that are neglected for simplicity in nearly all 
previous analyses of radiation damping.) 

Before proceeding with the analysis of (128!) . we would like to note for future use 
that the information we have gathered in the process of its derivation, and more 
specifically, equations ( l20l) and ( l25l) . allow Mikhailov's solution ( ITOl) to be rewritten 
purely in terms oi z = data as 

A first check on (l28l) is to note that it correctly reduces to md'^x^ / dr"^ = JF^ in 
the pointlike limit m ^ oo (where the Compton wavelength Zm ~^ 0). We will now 
perform a more substantial check by showing that it also makes firm contact with the 
results of [25] at finite m. 

In [25], two of us generalized the analysis of Mikhailov [30] to the case of a quark 
with finite mass, concentrating for simplicity on motion along one dimension. We 
showed that under such circumstances the total string (= field + quark) energy E 
and momentum P are no longer given by (fMj) and (fT6|) . but become 

with F the external force. These expressions show that for m < oo the rate (seen 
inside the integrals) at which energy/momentum is radiated by the quark differs from 
the Lienard result, and the formulas for the intrinsic energy Eq and momentum pq of 
the quark (given by the terms that follow the integrals) are similarly non-standard. 
Now, the total momentum P of the string (= field + quark) changes only due to the 
force that we exert on the endpoint (= quark), so dP/dt = F, or, using (l30!l . 

2.,n^[l-^F^)^dt[^^Z^.)- ■ ' ' 

Let us now compare this against our equation of motion (!28|) . For linear motion 
along direction x^, we have dx^/dr = ■j{l,v) and JF^ = ■y{Fv,F), so the = 1 
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component of (128|) reads 




(32) 



This is in precise agreement with (I3ip . Similarly, it is easy to see that the /i = 
component of fl28|) yields the expected result dE/dt = Fv with E as in fl30l) . 

We have thus verified that our quark equation of motion reproduces the en- 
ergy/momentum split between quark and radiation field previously deduced in [25] 
for the case of one-dimensional motion. It becomes clear then that ( l28l) encodes the 
covariant generalization of this split to the case of arbitrary motion — a result that 
would have been extraordinarily difficult to obtain using the non-covariant approach 
of [25]. To make this generalization explicit, we rewrite our equation in the form 

d dv^ d 

4-^-^ + -^ = ^'^ (33) 
dr dr dr 

recognizing 



as the intrisic four-momentum of the quark, and 



= dr 2^rn _ (^34^, 



dr 27rm2 \ 1 - ^-^-^^^ 



(35) 



as the rate at which four-momentum is carried away from the quark by chromo- 
electromagnetic radiation. 

Using again the fact that JF ■ d-rX = 0, we can immediately deduce from ( IMl) the 
mass-shell condition = — m^, which shows in particular that is indeed a four- 
vector, and so the split = Pg+Pj^^ defined in fl33l)- fl35l) is correctly Lorentz covari- 
ant. As we indicated above, P^^^ represents the portion of the total four-momentum 
stored at any given time in the purely radiative part of the gluonic field set up by 
the quark. The remainder, p^, includes the contribution of the near field sourced by 
our particle, or in quantum mechanical language, of the gluonic cloud surrounding 
the quark, which gives rise to the deformed dispersion relation seen in (134|) . In other 
words, is the four- momentum of the 'dressed' or 'composite' quark. All of this is 
completely analogous to the classical electromagnetic case that we briefiy reviewed in 
the Introduction, and in particular, to the covariant splitting of the Maxwell tensor 
achieved in [TT]. It is truly remarkable that the AdS/CFT correspondence grants us 
such direct access to this piece of strongly-coupled non-Abelian physics. 

Now that we have performed some checks on (l28l) and understood its proper 
physical interpretation, we should consider its implications. As noticed already in 
[25] for the case of linear motion, a salient feature of the equation of motion (l28l) . 
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as well as the dispersion relation (134|) and radiation rate (l35l) . is the presence of a 
divergence when JF^ = JF?^.;^, where 

•^crit = (36) 

is the critical value at which the force becomes strong enough to nucleate quark- 
antiquark pairs (or, in dual language, to create open strings) [24]. 

Let us now examine the behavior of a quark that is sufficiently heavy, or is forced 
sufficiently softly, that the condition ^ \\T'^\l2T\m^ < 1 (i.e., \T'^\ -C iJ^critl) holds. It 
is then natural to expand the equation of motion in a power series in this small param- 
eter. To zeroth order in this expansion, we have the pointlike result mci^x^/rfr^ = JF, 
as we had already mentioned above. If we instead keep terms up to first order, we 
find 

dr y dr 2T^m? j dr ^ ' 

In the 0(71) terms it is consistent, to this order, to replace with its zeroth order 
value, thereby obtaining 

(d^x^^ ^ d^x^\ _ ^IXd^x"" d^Xydx^ 
^ \d!?^ ~ 2^~d^ J ~ '2^~d^~d^~d^ ■ ^ ' 

Interestingly, ( 1381) coincides exactly with the Lorentz-Dirac equation ([2])! As ex- 
pected from the discussion in the preceding paragraphs, on the left-hand side we find 
the Schott term (associated with the near field of the quark) arising from the modified 
dispersion relation ( !34l) . On the right-hand side we see the radiation reaction force 
given by the covariant Lienard formula, as expected from the result ( fT^ [30] . which 
is the pointlike limit of the radiation rate ( l35l) . Moreover, by comparing ([2]) and (l38l) 
we learn that it is = V^/2vrm that plays the role of characteristic time/size tg 
for the composite quark. As we discussed around iQ, is the Compton wavelength 
(i.e., size of the gluonic cloud) of the quark, which is indeed the natural quantum 
scale of the problem. 

If we continue with the expansion of , the result at second order can be written 

as 

mx^ _ ^ «?M _ x^x^xn + (x"- 3F°f^x'^ - -rf^F ) = J^f" 



27r ' ' AiT^m \ 2 

We can of course continue this expansion procedure to arbitrarily high order in 
A/A|^^|/27rm^, and at order n in this parameter, we would obtain an equation with 
derivatives up to order -|- 2. Now, it is interesting to note that, in the case of 
non-relativistic QED, there are actually two different scales that appear in the cor- 
responding equation of motion [12]: while the terms with derivatives of fourth and 
higher order are all characterized by the quantum (Compton) radius Xc = h/m, the 
third-derivative (i.e., Abraham-Lorentz) term involves only the classical electron ra- 
dius cte = e^/mc^ -C Ac- In our strongly-coupled non-Abelian setting, the analogs 
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of these two scales happen to coincide. On the one hand, is analogous to Ac in 
that it gives the size of the cloud of virtual particles surrounding the quark, which as 
we explained at the end of Section [21 is set by the meson (and not the quark) mass, 
Zm = ^/tTi-^es- On the other hand, Zm is also analogous to cte in that it gives the radius 
of a charge distribution whose chromoelectrostatic energy equals the quark mass m, 
if we take into account the strong- coupling form of the potential V{L) oc \f\/ L |53j . 

As promised in the Introduction, our full equation fl2Sl) is thus recognized as an 
extension of the Lorentz-Dirac equation that automatically incorporates the size z^ 
of our non-classical, non-pointlike and non-Abelian source. The passage from (!38|) to 
(128|) . which can be viewed intuitively as the addition of an infinite number of higher 
derivative terms {z^d / dr)"^ , has a profound impact on the space of solutions. Here 
we will limit ourselves to two general observations, leaving the search for specific 
examples of solutions to the next section. The first is to notice that, unlike its 
classical electrodynamic counterparts ([1]) and ([2]), our composite quark equation of 
motion has no pre-accelerating or self-accelerating solutions. That is to say, the 
behavior of the quark at any given time r does not depend on jF^(r') at r' > r, and, 
in the (continuous) absence of an external force, (l28l) uniquely predicts that the four- 
acceleration of the quark must vanish. Our second observation, however, is that the 
converse to this last statement is not true: constant four-velocity does not uniquely 
imply a vanishing force. We will expand on this in the examples of the next section. 

Notice that, in the systematic approach from our result (!28l) to the Lorentz-Dirac 
equation, the pathologies that afflict the latter appear only in the very last step, when 
thex'^ term is introduced upon approximating (1371) by (l38l) . Indeed, it has often been 
advocated to eliminate these pathologies by a 'reduction of order' procedure (see, 
e.g., [SHE]), which treats the damping force as a perturbative correction and thereby 
justifies the replacement of fl38|) by fl371) . (A closely related hne of reasoning can be 
found in [55].) It is therefore satisfying to see that the dressed quark equation of 
motion predicted by AdS/CFT, Eq. ( l28l) . automatically comes out in 'reduced order' 
form. 

It is curious to note that (l28l) . which incorporates the effect of radiation damping 
on the quark, has been obtained from (ITOl) . which does not include such damping for 
the string itself. The supergravity fields set up by the string are of order and 
therefore subleading at large A''^. Even more curious [IB] is the fact that it is precisely 
these suppressed fields that encode the gluonic field profile generated by the quark, as 
has been explored in great detail (mostly at finite temperature) in recent years (see, 
e.g., [56] and references therein). It would be interesting to explore how the split into 
near and radiation fields is achieved from this perspective, but we leave this problem 
to future work [32] . 

It is natural to inquire whether the equation of motion (l28l) can be encoded in 
a variational principle. Since the complete system includes the composite quark in 
interaction with its radiation field, and the four-momentum of the latter is given by 
an integral over the quark worldline, at least naively we would expect that, if it is 
at all possible to write down an action, it should depend bilocally on the worldline. 
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This would rule out the obvious candidate action, namely S = S'ng + Sp with the 
classical solution (ITUi) plugged in. Indeed, the latter procedure simply leads again 
to (1181) (with r in the Nambu-Goto term relabeled to f), which (1231) allows to be 
converted into 

S = -m ! ^ + ! rfr A^(x(r))x^(r) (39) 

Unlike what happened in the pointlike (and no radiation damping) limit m — > oo, in 
this case variation of S with respect to holding J-'^ fixed does not yield the correct 
equation of motion fl28l) . Of course, S is by definition the correct on-shell action for 
the system, in the sense that it yields the right number when evaluated on a given 
worldline (using the J^^(t) obtained by solving fl2S]) ). but it would somehow need to 
be rewritten in bilocal form to explicitly show all relevant dependence and thus 
constitute the desired variational principle. 



3 Examples 

3.1 One-dimensional motion 

We have already seen in Section that, in the case of one-dimensional motion, our 
general equation of motion (1251) reduces to fl52]) . The latter can be further simplified 
to 

where we have preferred to express the prefactors in terms of the quark Compton 
wavelength instead of its mass m, and as in fl25l) we have used the abbreviation 
f = (27r/v^)F. Here we see directly that, in contrast with the usual case, the 
acceleration at any given time depends not only on the value of the applied force but 
also on its rate of change. 

It follows from (l40l) that a quark that is free for any extended period of time will 
not accelerate. In other words, as we had already indicated in the general discussion 
below (1551) . there are no self-accelerating solutions, which is just as one would have 
expected given the extended nature of the charge. On the other hand, the very fact 
that the charge has a 'deformable' internal structure implies that there is more than 
one way to get it to follow any given trajectory. Indeed, for any choice of f (t), POl) 
fixes the applied force ^{t) not algebraically, but through a differential equation that 
inevitably gives rise to a one-parameter family of solutions (differing by their initial 
conditions). 

The simplest example of this non-uniqueness of the external force is the case of 
constant velocity, where (l40l) can be easily seen to imply that 

f{t) = ±\^ech(^^^^{t-t,)\ , (41) 
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with to an integration constant. Only for ~^ ioo does one recover the simple 
result F{t) = 0. For all finite values of the integration constant, the force starts out 
at zero at asymptotically early times, rises steadily until it attains the critical value 
given by fl36l) a.t t = to, and then approaches zero again as t — > oo. 
It is certainly peculiar that one can continually apply a force to the quark and still 
have it move at constant velocity. Nonetheless, it is easy to verify through numerical 
integration that indeed the application of the force (HTj) to the endpoint of a string 
whose initial profile is chosen in accord with fl2^ produces no acceleration. The 
energy provided to the system by F{t) does not translate into an increase of the 
endpoint velocity, but into a continuous modification of the string tail, or, in gauge 
theory language, a change of the gluonic field profile. In fact, with the formulas 
derived in the previous section, we can make a much more precise statement: ([3] 
and fl35l) reduce to [251 




7m 



rad 



dt 



where we see that application of the force fHT]) results in 



dE^^^ 
dt 



{coth ( ^ ] =p i; csch ( ^ ] 



7m 



csch 



t-k 



1 — V sech 



t - tc 



which show precisely what fraction of the energy goes into (or comes out from) re- 
arranging the near gluonic field, and what fraction goes into radiation. In the limit 
to ±00, no force is applied and we of course recover Eq = 7m, dE^^/dt = at all 
times. Similar conclusions can be drawn about the linear quark momentum. 

Notice that the width of the time interval over which the force f HT]) differs appre- 
ciably from zero is just the Compton wavelength of the quark, Zm (with an appropriate 
Lorentz dilation factor), which serves to emphasize that this peculiar phenomenon is 
made possible only due to the non-pointlike nature of our non-Abelian source. 

This non-uniqueness of the force is equally present for arbitrary trajectories. An- 
other example that is easy to study is the case of constant force. If we assume that 
Fit) = constant, (l40l) reduces to d{^v)/dt = Zm^/ \/l — z^f^, which has the same 
form as the standard equation of motion for a particle with constant proper acceler- 
ation, except that the force here is non-linearly rescaled. The solution is thus 



x{t) = Xq ± 



1 



-zl+it- tof 



where Xq and to are integration constants. But if we run the argument in re- 
verse, and ask what type of force would lead to the hyperbolic motion x(t) = 
Xq lb -y/C^ + (t — to)^, we see again that the differential equation (l40l) admits solutions 
other than the obvious ^(t) = 1/zm^/C^ 
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Before closing this subsection, we would like to make one additional observation. 
After some straightforward algebra, it is easy to see that (1351) can be rewritten in the 
form 



"-'^rad 
dT 



^2 



1 - 



A 



jr2 



(42) 



Plugging this into 0331) . one obtains a first order differential equation for with the 
same structure as the Langevin equation, but with a friction coefficient that depends 
on the external force, 

dp^ 
dr 

where /i(J^) = {y/X/2Trm^)J^y ^1 - X/An^m^T"^. In the case of motion in one spatial 
dimension, this equation can be solved analytically. 



(43) 



Pq = exp 



-dx 



A + 



exp 



27rm^ 



1 - 



A 



-dx ) T{y)dy 



(44) 



with A an integration constant. We should stress that this solution is valid only 
in the case of one dimensional motion, where JF^ = F"^. In two or three spatial 
dimensions we have not been able to construct an explicit solution, because in that 
case JF^ involves the velocity of the quark, and remains undetermined until we find the 
quark trajectory. Nevertheless, it is interesting to note that, in the case of constant 
JF^ (which is not the same as constant F^), equation (l43l) can be interpreted as a 
Langevin equation of motion for a particle with momentum pq (up to a constant 
term). The force in terms of the physical time t is just the force needed to move the 
quark in a 'dissipative medium' characterized by the constant friction coefficient /x. 

It would be interesting to explore the dynamics of the dressed quark in two or three 
spatial dimensions. However, as we have just remarked, in that case fl28|l is highly 
nonlinear not only in the external force (which it was already in one dimension), but 
also in the velocity of the quark. For this reason, it is unfortunately very difficult to 
find analytic solutions to it. In the heavy quark (or small force) approximation where 
(128|) linearizes and reduces to the Lorentz-Dirac equation (!38l) . we could of course 
carry over to our setting the various solutions that have been worked out in the 
past. For instance, in the first reference of [7j, Rohrlich was able to find an analytic 
solution for a central force problem using the Frenet equations. Starting from this 
and perturbing the solution it should be possible to obtain (at least numerically) 
the first correction beyond Lorentz-Dirac predicted by our framework, and deduce 
for instance the r ate of synchrotron radiation. Such a calculation might shed some 
additional light on the physics behind these extended objects, but is not central for 
the purposes of our analysis here, so we prefer to leave it for future work. 
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3.2 Nonrelativistic limit 



Let us now choose a specific Lorentz frame {x,t), and restrict attention to motions 
such that the quark velocity v = dx/dt, acceleration a = d'^x/dt'^, jerk J= d^x/dt^, 
and all higher derivatives, as well as the force ^ and its rate of change d^/dt, are 
small in units of the Compton wavelength Zm- Under such conditions dr ^ dt, and 
the spatial components of TI7}i adopt the linearized form 

d oc d X '-yd 

d^-'-^ + ^m^---- = 'rnr- (45) 

Adding to this expression its time derivative multiplied by Zm, we obtain the simplified 
form 

d'^x , dF . . 

^ = + (46) 

which is the linearized version of (128|) . We see here very directly that, as explained 
in Section \2.2\ the unusual dependence on the rate of change of the force encodes an 
infinite number of higher derivatives of the quark trajectory x(t), derivatives which 
in turn refiect the extended nature of the quark. 

It is also useful to note that in this linearized limit, the string embedding fl29l) can 
be rewritten in the form 

X{t, Z) = f(tret) + Zm) X] ("^m)'"^ ^ (^ret) , ^ret = t - Z + Zm , (47) 

1=1 

which involves only the quark worldline x{t) and not the force F(t) . If we plug this into 

the (quadratic version of the) total string energy E{t) = — J dzH^ = J dz{X'^+X''^)/2, 
and imitate Mikhailov's procedure [20] (see Section I2.2p . the expression naturally 
splits into 

E{t)= f dt,^,^ + E,it) , 



with 




and 

/ oo \ ^ / oo 

d Ej i^Qj^ I ^ ^ , , ^ 2 d X , , \ ^TTi dx , , / ^ ^ , > ^ 2 d X 



Using dlSD, it is easy to check that these expressions indeed coincide with the non- 
relativistic limit of the time component of fl5^ and fl35l) . The spatial components 
can be verified in a similar manner. Notice that, in the relativistic case analyzed 
in [25] and Section 12.21 of the present work, the explicit presence of the force in 
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the solution (l29ll prevents us from deriving the spht (133|) as we did here, directly 
imitating Mikhailov's procedure. It is therefore comforting to see that, at least in 
the non-relativistic case, the direct result a.t z = indeed agrees with what we had 
previously inferred indirectly from Mikhailov's unambiguous split for the auxiliary 

data at z = 0. 

Just like in the fully relativistic setting, if we directly substitute P7j) into the 
(quadratic version of the) string action -\- we do not arrive at a variational 
principle that correctly encodes the equation of motion fHS]) or fHBl) . Here, however, 
it is easy to write down an action that does the right job. In fact, we have not one 
but (at least) two options: we can evidently get (HSj) from 



Sr. 



m 



dt 



dx\ ^ / d'^x\ 



d^x\ 



+ dtF-x 



and (H^ evidently follows from 



dt'm-^im 



It is curious to note that the action S^r in terms of higher order derivatives of 
the position vector with respect to r is reminiscent of the dynamical description of 
a particle in noncommutative symplectic mechanics given in [57] , where it was found 
that the action can be written in terms of higher order derivatives of the position 
vector despite the fact that the equations of motion are of second order. 

Let us now explore some solutions. Either from S'^^ or directly from the equation 
of motion (H6l) . one immediately has a first integral of the motion. 



dx 
'di 



dt'm+4nm 



The case with zero acceleration corresponds to f + z^,^ = 0, whose solution is 
fit) = f^Qexp[— (t — to) /-2m.]- This is of course simply the non-relativistic limit of 
(141 p . Since we are working now in the linearized approximation, this solution of the 
homogeneous equation of motion can be added to any solution of the inhomogeneous 
equation fH6|l to yield another solution, so it is very clear that the non-uniqueness 
of the force is present for any quark trajectory whatsoever. The reason is by now 
familiar to us: the energy provided to the dressed quark by this one-parameter family 
of forces has the effect of modifying the gluonic field profile, and consequently does 
not translate into an increase of velocity. 
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